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DIFFERENTIAL ANALYSIS FOR THE TURBULENT BOUNDARY LAYER ON A 
COMPRESSOR BLADE ELEMENT (INCLUDING 
BOUNDARY-LAYER SEPARATION) 
by James F. Schmidt and Carroll A. Todd 
Lewis Research Center 

SUMMARY 

In order to provide an approximate boundary- layer prediction at the trailing edge of 
a compressor blade (especially with separated flow), a two-dimensional differential 
analysis is developed for the turbulent boundary layer (including the separated reverse- 
flow region). Prandtl's classical mixing- length theory is used with a single eddy- 
diffusivity model throughout the boundary- layer region. The use of a new shear- stress 
function (instead of a velocity function) for the dependent variable in the momentum equa- 
tion makes it possible to calculate mathematically stable solutions into the boundary- 
layer separation region. 

This analysis is applied to a blade in a cascade with three different inlet-flow angles 
for comparison with experimental boundary- layer measurements. The predicted 
turbulent- boundary- layer thicknesses and velocity profiles are in good agreement with 
experimental data, even in the separated region. For the thickest turbulent boundary 
layer on the blade with a separated region, 8 minutes of computing time (Lewis DCS 
7094-7044) are necessary for this calculation. 


INTRODUCTION 

In the analysis of compressor blade performance, many basic flow parameters such 
as maximum blade loading, blade loss, and blade stall are primarily related to the 
blade-surface boundary layer. Therefore, the need exists for a boundary- layer calcula- 
tion procedure to aid in the blade- row flow analysis, especially in the flow- separation 
region. To select or develop a suitable boundary- layer calculation technique, some 
characteristics of the flow through compressor blade rows must be known. 



All compressor blade rows are designed to increase the static pressure of the fluid. 
Thus, the blade- surface boundary layer experiences an adverse pressure gradient. 

With high blade loading and near-stall operation, these adverse pressure gradients be- 
come very strong and a separated boundary- layer flow region often occurs. 

In general, the airflow entering the blade rows is very turbulent (high turbulence 
intensity level) because it has passed through an upstream blade row and a long inlet. 

The shear wakes from the upstream blade row and the strong adverse-pressure-gradient 
flow history produce a flow environment of high free-stream turbulence intensity. 

Over much of their operating ranges, compressor blades produce large velocity 
peaks very close to the nose of the blade (on the suction surface). And thereby, adverse 
pressure gradients are expected over almost the entire length of the blade. These 
strong adverse pressure gradients (shortly after the leading edge of the blade), as well 
as high turbulence intensity levels, tend to cause early transition or separation and 
reattachment as a turbulent boundary layer. Therefore, the turbulent boundary layer 
effectively exists over the entire blade. Hence, a theoretical analysis is needed which 
can calculate the turbulent boundary layer (with high turbulence intensity levels) through 
the blade stall flow- separation region. However, a turbulent-boundary- layer calcula- 
tion in the separated flow region is not only mathematically very difficult but also de- 
pends upon the unknown and unmeasured Reynolds turbulent shear stress. Thus, this 
boundary- layer prediction must be considered an approximation. 

In the past, either a turbulent-boundary- layer growth correlation or, more recently, 
an integral boundary- layer analysis (ref. 1) was used to predict the boundary -layer 
growth on a compressor blade element. Most integral analyses, however, do not have 
the capability of predicting the boundary- layer growth in the separated reverse-flow 
region. 

In addition, several differential numerical calculation methods have been developed 
for the turbulent boundary layer (refs. 2 to 6). All these calculation methods use an 
eddy-diffusivity concept for the unknown turbulent Reynolds shear stress and a finite- 
difference calculation procedure which utilizes a velocity function as the prime depend- 
ent variable across the boundary layer. Considering that the velocity profile will have 
an inflection point at and beyond the separation point, the authors believe the use of a 
velocity function as the prime dependent variable may cause large numerical errors and 
instabilities in the basic finite- difference differentiation technique. 

Presently, very few analytical methods are available for calculating turbulent 
separated flows caused by adverse pressure gradients. One such method (ref. 7) pre- 
sents an integral turbulent- boundary- layer analysis which was applied to turbulent flows 
near and in the separated region. However, a comparison of the integral analysis of 
reference 7 with experimental data (near and in the separated region) still shows con- 
siderable discrepancies in the boundary- layer thicknesses (ref. 7). In all these 
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turbulent-boundary- layer analyses (refs. 1 to 7), the free- stream turbulence level is 
assumed to be small; and this, in general, is not indicative of the boundary- layer flow 
over compressor blades. 

This report develops an analysis for the turbulent boundary layer on compressor 
blades, including the separated region with reverse flow. The present two-dimensional 
analysis is a successive- approximation, nonsimilar calculation (as in ref. 8) of the 
turbulent, differential, boundary- layer equations. The calculation system of the present 
analysis is suitable for calculating the turbulent boundary layer in strong adverse- 
pressure-gradient regions and also in reverse-flow separated regions. The use of a 
new shear- stress function for the dependent variable in the momentum equation makes 
it possible to calculate mathematically stable solutions into the separated region. 

At the separation point and in the reverse-flow separated region, even the complete 
boundary- layer equations (ref. 9) may not be valid. In addition, the eddy diffusivity is 
presently unknown in the separated boundary- layer region. Also, the momentum equa- 
tion normal to the blade surface is not included in the present system of boundary- layer 
equations, and the pressure is assumed to be constant across the boundary layer. For 
these reasons, the present turbulent- boundary- layer solutions in the separated flow 
region should be considered an engineering approximation. 

The present analysis also uses the eddy- diffusivity concept for the unknown turbu- 
lent Reynolds shear stress, but with a different eddy- diffusivity model than those used 
in previous analyses. Since high free-stream turbulence levels as well as extreme, 
adverse pressure gradients are characteristic of flows over compressor blades, a large 
free-stream shear flow is believed to describe the rotational inviscid flow over com- 
pressor blades. Thereby, the large shear stress in the turbulent boundary layer does 
not diminish to the usual insignificant value at the edge of the boundary layer. In fact, 
the shear stress at the edge of the blade turbulent boundary layer is assumed in this 
analysis to be a large significant fraction of the maximum shear stress in the boundary 
layer. This large- shear- stress behavior at the edge of the boundary layer has experi- 
mentally been measured in a turbulent boundary layer with high levels of free-stream 
turbulence intensity (ref. 10). Based on these reasons, Prandtl’s mixing- length theory 
(ref. 9) for the eddy diffusivity, with Van Driest's damping factor (ref. 11), is used 
throughout the entire boundary layer. 

The results of this analysis are compared with the experimental turbulent-boundary- 
layer measurements of reference 12. Reference 12 gives two- and three-dimensional, 
turbulent-boundary- layer measurements on a blade in a cascade for several inlet-flow 
angles. These measurements include the turbulent-boundary- layer, flow- separation 
region. The turbulence intensity level for this cascade flow is reported to be high 
(ref. 12). 
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ANALYSIS 


The present analysis is a nonsimilar differential calculation of the turbulent- 
boundary-layer equations. 

Assumptions 

The analysis is based on the following assumptions: 

(1) Surface curvature is neglected in the boundary- layer equations. 

(2) Flow is two dimensional and steady. 

(3) Normal Reynolds stresses are negligible, even for reverse flow. 

(4) Static pressure is constant across the boundary layer. 

(5) The eddy-diffusivity concept can be used for the tangential Reynolds shear 
stress. 

(6) Prandtl’s eddy-diffusivity, mixing- length theory is assumed to be applicable for 
the turbulent boundary layer on blade surfaces with very strong, adverse pressure gra- 
dients and high turbulence levels. 

(7) The eddy diffusivity is always positive, even for reverse flow. 

(8) The turbulent boundary layer begins at the stagnation point on the blade surface. 

Basic Equations 

In reference 8, the turbulent-boundary- layer equations for axisymmetric steady 
flow are presented in terms of the eddy-diffusivity concept. For two-dimensional flow, 
the system of boundary- layer equations in reference 8 reduces to 

Continuity: 

(PU) S + (PV) S - 0 (1) 

Momentum : 

PUU S ♦ pVUy = P e U e (U e ) g ♦ [<P ♦ Pduj^ (2) 


4 



Energy: 


PUHq + pVH - 

(— + - £ ^-\ H v 

+ 

& y 

V Pr L P-t/ y _ 




State: 


P = ptftt 


(3) 


(4) 


(All symbols are defined in appendix A. ) 

The boundary- layer coordinate system along a blade cross section is shown in fig- 
ure 1, where S is the boundary- layer coordinate parallel to the blade surface, c is 
the chord length, x is the distance coordinate from the nose to the trailing edge of the 
airfoil section, and x/c is the percentage-of- chord location (0 to 100). 


Eddy-Diffusivity Model 


Before the turbulent-boundary-layer equations (1) to (4) can be solved, an expres- 
sion for the eddy diffusivity e is needed. As discussed in the INTRODUCTION, a one- 
layer model is assumed for the eddy diffusivity. The eddy diffusivity is given by Van 
Driest's modification (ref. 11) (viscous damping effect near the wall) of Prandtl's 
mixing- length expression (ref. 2): 


e 


2 2 
= k y 




(5) 


where A is (from ref. 2) 


A 



( 6 ) 


and k is an experimental flow constant (0. 40) based on low-speed incompressible data. 
The dynamic viscosity of air is approximated by 
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Transformed Equations 


Using a compressibility type of transformation (ref. 13) yields the following trans- 
formation equations, used in solving the boundary-layer equations (1) and (2): 


P e (U e ) 1/2 r P 
V =—5—5 / — dy 

«Vo 2S > 1/2 7o Pe 

(8a) 

s = s 

(8b) 

The partial derivatives transform as follows : 


P U g /2 

(, y= - >I/2 ( V () s s y V° 

(p 0 M 0 2S) I/z 

(8c) 

( ) s = ( >s + ( Vs 

(8d) 

From continuity, the stream function \Jy is defined such that 


pu = (^) y 

(8e) 

- pv = -(ijs) s 

(8f) 

By letting 


* =(p 0 !* 0 2SU e ) 1/2 f(rj,S) 

(8g) 

Then, 


p| 

ii 

(8h) 


Applying the transformation equations (8) to the momentum equation (2) and reducing 
give the following form for the momentum equation: 
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C + Mt) 2 £ V + ( ®2 + = 2Sf A) S ' 13 -f 

\ P e/ J Jri — 


- (f Y 


where the parameters in equation (9) are as follows: the density- viscosity ratio C is 
defined as 


P T 

pp e Op 


P 0 P 0 P 0 t M 0 


Using equations (7) and (8f) with the definition of total temperature 


T = t + 


in equation (10) and reducing result in the transformed density-viscosity ratio 




p o i - D 3 (y 2 


and the ratio of kinetic to total energy 


2c p T 0 


The density ratio p/p e can be expressed as 


„ t 1 - D. 
p _ e 3_ 


p e 1 


1 - °3 V 
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The flow coefficients and D2 are 


p e 

Dj = — - — (14) 

P 0 P 0 

D, = 1 + ^ (15) 

1 2 


where the velocity gradient parameter j3 is 


J8 



(16) 


And S is the boundary- layer distance coordinate starting at the origin of the turbulent 
boundary layer. 

In the momentum equation (9), the terms involving the velocity derivatives with 
respect to S, 2Sfgf^, and 2Sf^(f^)_ are the nonsimilar terms that were neglected in 


the similar- solution calculation of reference 8. The energy equation is transformed and 
programmed in the computer program but is not indicated in the present analysis. This 
omission is believed to result in a less complex presentation of the theory, especially 
since the test cases do not include energy transfer. 

With substitution of the transformation equations (8) into equations (5) to (7) and 
reduction, the transformed eddy diffusivity becomes 




^e/£f^V /2 

"w V 2 P J 


(17) 


( 18 ) 
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Grouping the molecular viscosity and eddy diffusivity into one coefficient gives 


C 


1 


C + Dj 



£ 


( 19 ) 


Substituting equation (19) into equation (9) yields the following momentum equation: 


( C lO + <fD 2 + 2Sf S )f „„ = 2§ W5 - P 


P 


(f Y 


( 20 ) 


Multiplying the momentum equation (20) by the velocity gradient parameter j3, which is 
constant across the boundary layer, gives 


K l rrn) * < fD 2 * 2Sf S^V = ^ - P 
T 1 ^ 


Now the momentum equation (21) has the proper form such that the following new depend- 
ent variable transformation can be used: 

& = flCj f^ (22) 

This substitution (eq. (22)) eliminates the need for numerically differentiating the eddy 
diffusivity. And more importantly, this transformation is believed to be an important 
factor in stabilizing the numerical boundary- layer solutions, especially for the adverse 
pressure gradients approaching flow separation and beyond. 

When equation (22) is substituted into equation (21), the momentum equation in 
terms of the new dependent variable YP has the form of an ordinary differential equation 

YP + a^ = h (23) 


- - ( f n> 2 
P V 


( 21 ) 


where 
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a = 


(24) 


fD 2 + 2Sfg 


h = * 2S ^Vc - ^ 


— - (f ) 2 

_P V 

L P e 


(25) 


The following boundary conditions are applied to the momentum equation (23): 


f = f = 0 at rj = Cp 
V 

f 77 = 1* 0 at * = ^ej 


(26) 


The nonsimilar streamwise velocity derivatives (fg and (f )_) in equations (24) and 

(25) are evaluated with a simple two-point numerical derivative formulation. This 
simple formulation of the streamwise velocity derivatives reduces the momentum equa- 
tion (23) to a linear, ordinary, differential equation (function of 17 only) and is easily 
integrated in closed form. The detailed numerical calculation method is given in appen- 
dix B. 


Skin- Friction Coefficient and Boundary- Layer Thicknesses 

The following definitions of various boundary- layer terms are useful for analyzing 
the behavior of the boundary layer, as well as for comparison with experimental data: 
The shear-stress distribution across the boundary layer is expressed as 

t = (M + Pe)-— (27) 

3y 

When equations (8), (19), and (22) are substituted into equation (27), the nondimensional 
shear-stress function becomes 


£ 

0 


<PqMo 

I 2S , 


1/2 


U 


3/2 


( 28 ) 
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From equation (28) the shear-stress function # is shown to be the nondimensional 
shear-stress distribution across the boundary layer. At the wall, the shear stress re- 
duces to 


r = u (U = - p U 2 C* (2 

w 'wV y/ n e f 

Substituting the transformation equations (8) into equation (29) and solving for the skin- 
friction coefficient C f give 



/vr 

\U 2S / 

\ e / 



(30) 


The following boundary- layer thicknesses are very useful for comparison with ex- 
perimental measurements and other theoretical analyses: 



TURBULENT REYNOLDS SHEAR-STRESS MODEL 

Before a comparison of the present theory with experimental data is made, a de- 
scription of the present turbulent Reynolds shear- stress model and how it differs from 
other theoretical analyses is helpful. Most theoretical analyses use the eddy-diffusivity 
concept which states that the turbulent Reynolds shear stress is equal to the eddy dif- 
fusivity times the velocity gradient across the boundary layer. 

Several analyses such as references 2, 3, and 6 use Prandtl's mixing- length theory 
with viscous damping near the wall matched to an outer constant eddy-diffusivity (wake) 
region. A slightly different eddy-diffusivity model, as in reference 5, consists of 
Prandtl’s mixing- length eddy diffusivity for the entire boundary layer (like eq. (5)) but 
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with two different mixing- length regions. The first mixing- length region of reference 5 
is Prandtl's classic mixing- length expression (Z = /cy), and the second region is a con- 
stant approximated from experimental data, as shown in figure 2. The measured, non- 
dimensioned, mixing length Z/5 greatly increases near the edge of the boundary layer 
for low free-stream turbulence levels in adverse-pressure-gradient flows (refs. 14 and 
15). However, the mixing length in reference 5 was assumed to be a constant in this 
region of the boundary layer. 

Figure 3 gives a generally good comparison of Prandtl's mixing length with the 
general trend of the measured mixing length (ref. 9) for flows with high free-stream 
turbulence levels. Because of this favorable comparison of Prandtl mixing- length 
theory with data and the fact that flows with high free-stream turbulence levels generally 
occur over compressor blades, Prandtl's mixing- length theory is used for the eddy- 
diffusivity model in the present analysis. 


COMPARISON OF PRESENT THEORY WITH EXPERIMENTAL MEASUREMENTS 

The present turbulent-boundary- layer analysis was compared with the experimental 
boundary- layer data of reference 12 for a straight cascade at several inlet-flow angles. 
The cascade of reference 12 consisted of nine blades with a chord length of 12. 4 centi- 
meters, a span of 50. 8 centimeters, and an NACA 65-410 airfoil section. The free- 
stream flow properties were calculated from the static- pressure measurements along 
the airfoil that are presented in reference 12. The boundary- layer measurements on the 
cascade airfoil were compared with the present analysis for the following three different 
inlet- flow conditions: 

(1) Two-dimensional flow with a cascade inlet-flow angle of 55° (2D-55 0 ) 

(2) Two-dimensional flow with a cascade inlet-flow angle of 65° (2D-65 0 ) 

(3) Three-dimensional flow with a cascade inlet-flow angle of 60° (3D- 60°) 

In the third flow condition (3D- 60°), three-dimensional flow was simulated by pro- 
ducing an added velocity diffusion near one end of the cascade blades. (Although the 
present analysis is strictly a two-dimensional analysis, the present calculation system 
was still applied to this simulated three-dimensional flow. ) The added velocity diffusion 
was artificially produced by placing an obstruction downstream of the span-wise center- 
line of the cascade (ref. 11). The turbulence intensity level for this cascade is reported 
to be high (ref. 12). For all inlet- flow conditions, the flow over the cascade blades is 
incompressible. 
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Boundary- Layer Thicknesses 


In figures 4 and 5 the predicted momentum and displacement thicknesses along the 
airfoil section are compared with experimental data (ref. 12) for the three cascade 
inlet-flow conditions. In general, the predicted momentum and displacement thickness 
distributions along the airfoil are in favorable agreement with the experimental data. 

For the suction surface with the inlet-flow condition 2D- 55°, the greatest discrep- 
ancy (10 to 20 percent) between the predicted and measured thicknesses occurred at 80 
and 90 percent of chord (fig. 4(a)). This discrepancy is believed to result partially 
from an inaccurate velocity gradient /3 distribution along the airfoil. Since the 
boundary- layer growth is a strong function of |3, any error in 0 will significantly affect 
the boundary- layer thicknesses. The calculation of j3 involves a numerical differentia- 
tion of an experimentally determined free -stream velocity distribution. Thereby, two 
sources of error are apparent: the reduction of experimental pressure measurements 
to a velocity distribution and the numerical differentiation of this velocity distribution. 

In addition, reference 12 presents the pressure measurements in graphical form with- 
out a grid, which necessitates an interpolation of measured results - another source of 
error. 

Figure 4(b), for the pressure surface 2D-55 0 , shows that the predicted thicknesses 
reduced more rapidly than measured values as the favorable pressure gradient was en- 
countered at the last 25 percent of the chord. Again, this difference between the pre- 
dicted and measured thicknesses is believed to be caused primarily by an inaccurate 
distribution of |3 along the airfoil. 

Although the turbulent boundary layer separated on the suction surface for inlet-flow 
conditions of 2D- 65° and 3D- 60°, the predicted thicknesses still compared favorably 
with measured values (figs. 5(a) and (b)). The distribution of the predicted skin- friction 
coefficient along the suction surface is also given in figures 5(a) and (b). The point on 
the airfoil surface at which the skin- friction coefficient first becomes negative is a good 
indication of the location of boundary- layer separation. Figures 5(a) and (b) show that 
the turbulent boundary layer separated at about 70 percent of chord for the inlet- flow 
condition 2D- 65° and at approximately 40 percent of chord for the inlet-flow condition 
3D- 60°. 


Velocity Profiles 

In figures 6 to 9 the predicted velocity profiles are compared with experimental 
data for the three cascade inlet-flow conditions at 20 or 98 percent of chord. These 
velocity profiles are presented as a function of the transformation variable tj, which is 
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independent of the skin-friction coefficient. Usually, the velocity profile is presented 
as a function of a wall distance variable y + (ref. 16), which requires a knowledge of 
the skin-friction coefficient. Since t? does not depend upon the skin-friction coefficient, 
the present comparison of predicted velocity profiles with experiment data does not con- 
tain any error from the experimental or predicted skin-friction coefficients. Also, near 
the separated region the skin-friction coefficient becomes very small and experimentally 
susceptible to large errors (if determinable at all). 

For the first cascade inlet-flow condition (2D-55 0 ), the predicted velocity profile at 
20 percent of chord is in excellent agreement with the limited experimental data (fig. 6). 
Even at 98 percent of chord, the predicted velocity profile is still in good agreement 
with experimental data (fig. 7) for the same inlet- flow condition (2D- 55°). For a higher 
cascade inlet-flow angle (2D- 65°) with increased velocity diffusion, figure 8 shows good 
agreement of the predicted velocity profile with experimental data at 98 percent of chord. 
Similarly, for even larger velocity diffusion (3D- 60°), the predicted velocity profile is 
still in good agreement with experimental data at 98 percent of chord (fig. 9). This 
good comparison of the predicted with the experimental velocity profiles shown in fig- 
ures 8 and 9 is of special interest because the turbulent boundary layer had already 
separated upstream of the 98-percent-of-chord station. In fact, for the largest velocity 
diffusion case (3D-60°), the predicted turbulent- boundary- layer separation point 
occurred at approximately 40 percent of chord. Therefore, the present turbulent- 
boundary-layer theory is shown to predict accurate velocity profiles and thereby accu- 
rate boundary- layer thicknesses, even for a significantly separated turbulent boundary 
layer on this particular cascade blade. However, the present theory has thus far only 
been applied to this one set of cascade data for a separated, turbulent, boundary layer. 
Additional separated- boundary- layer dgta over a wider range of flow conditions are 
needed to assess the accuracy, generality, and applicability to different flow orienta- 
tions of the present analysis. 


CONCLUSIONS 

In order to provide an approximate boundary- layer prediction at the trailing edge of 
a compressor blade (especially with separated flow), a two-dimensional differential 
analysis was developed for the turbulent boundary layer (including the separated reverse- 
flow region). The use of a new shear-stress function (instead of a velocity function) for 
the dependent variable in the momentum equation made it possible to calculate mathe- 
matically stable solutions into the boundary- layer separation region. Prandtl's classical 
mixing-length theory was used with a single eddy-diffusivity model throughout the 
boundary- layer region. 
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This analysis was applied to a blade in a cascade with three different inlet-flow 
angles for comparison with existing experimental boundary- layer measurements. The 
predicted turbulent-boundary- layer thicknesses and velocity profiles are in good agree- 
ment with experimental data, even in the separated region. For the thickest turbulent 
boundary layer on the blade with a separated region, 8 minutes of computing time 
(Lewis DCS 7094-7044) are necessary for this calculation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 20, 1973, 

501-24. 
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APPENDIX A 


SYMBOLS 


A constant, eq. (6) 

a term in momentum equation, eq. (24) 

C density-viscosity ratio, eq. (11) 

Cj skin-friction coefficient, eq. (30) 

transformation parameter, eq. (19) 
c chord length 

Cp specific heat at constant pressure 

flow coefficient, eq. (14) 

D2 flow coefficient, eq. (15) 

Dg flow coefficient, eq. (12) 

dU e 

velocity derivative at edge of boundary layer 

dS 

e natural or Naperian base (e = 2. 71828) 
f velocity ratio, U/U e 

H total enthalpy 

h term in momentum equation, eq. (25) 

Iq iteration criterion, eq. (B16) 

Kj constant of integration, eq. (B15) 

K 3 constant of integration (Kg = 0) 

Kg constant of integration (Kg = 0) 

l mixing length 

P static pressure 

dependent variable transformation, eq. (22) 

Pr k laminar Prandtl number 

Pr^, turbulent Prandtl number 

<R universal gas constant 

S boundary- layer coordinate parallel to blade surface, S = S 
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T total temperature 

t static temperature 

U velocity parallel to blade surface 

V velocity normal to blade surface 

x distance coordinate from nose to trailing edge of airfoil section 

y boundary- layer coordinate normal to blade surface 
y + wall distance parameter, eq. (18) 

|3 velocity gradient parameter, eq. (16) 

5 boundary- layer thickness 

5* displacement thickness 

e eddy diffusivity 

r\ transformed coordinate in y-direction, eq. (8a) 

6 momentum thickness 

k empirical constant (0. 40) 

M molecular viscosity 

p density 

r shear stress 

\f/ stream function 

w integration factor, eq. (Bll) 

Subscripts: 

e edge of boundary- layer condition 

Q iteration number (increase in T) e ) 
w wall condition 

y, yy derivatives of function with respect to y‘ or S coordinate 

77,7717 derivatives of function with respect to the transformed rj coordinate 

0 reference stagnation condition 

Superscript: 

d exponent for viscosity variation with temperature (0. 65 for air) 



APPENDIX B 


CALCULATION PROCEDURE 


For ease of reference, the transformed momentum-boundary- layer equation and 
the boundary conditions are repeated here. 


+ a c ^ = h 
V 

where 


fD 2 + 2Sfg 



(Bl) 


(B2) 


h = e 2S \(% - P 


- (f r 


(B3) 


with the boundary conditions 

f = f = 0 at r? = 0 (B4) 

\ • 

f^ = 1.0 at tj = rj e (B5) 

where rj is the boundary- layer thickness in the transformed tj- coordinate. 

The nonsimilar streamwise derivatives fg and (f^)_ are approximated in the fol- 
lowing manner: 



f " %-l 


AS 



( f ri)g - ( f r?)s_! 


AS 


(B6) 


(B7) 
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where fg j and (f^)g j are the velocity functions at the previous S location. The 
nonsimilar derivatives are underrelaxed as follows: 






(B8) 




L MJR 


= (f ) + 0. 1 

Kr »s-1 


( f T))g " Wg-i 


(B9) 


where the subscript R denotes relaxed derivatives. 


Solution of Momentum Equation 

The simple two-point formulation of the streamwise velocity derivatives, equa- 
tions (B6) and (B7), reduce the momentum equation (Bl) to a linear, ordinary, differen- 
tial equation (function of p only). This linear, ordinary, differential equation (Bl) of 
the first order can readily be solved in the following manner: 



cuh dp + 


K-, 


w 


where 


and 




(BIO) 


(Bll) 


(B12) 


(B13) 


(B14) 
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When the boundary conditions (B4) and (B5) are applied to equations (BIO) to (B14), 
the constants of integration (Kj, Kg, Kg) become 



The general procedure for determining the edge of the boundary layer (correct r? g ) 
is to compare the velocity derivative at the wall (f ) between successive increments 


w 


in T 7 g . If the iteration criterion for the change in (f ) is met, the edge of the bound- 
ary layer is obtained. In detail, the correct tj 6 is determined by the following iteration 
procedure: 

(1 ) Assume some initial tj and increase « by An for each outer iteration. 

6 c 6 

(2) Use the following iteration criterion: Let 





Q+i L 


(W 


w 


Q 



(B16) 


where Q indicates the iteration number (increase in rj e ). If |Iq| < 0. 005, the correct 
tj and boundary- layer solution are obtained. 


General Flow of Program 
The general flow of the program is as follows: 

(1) Input the total flow conditions (Pq, Tq, Pq, etc. ) and the table of S against 

V P 0' 

(2) For a given distribution of S (interval in S), calculate all free-stream condi- 
tions, such as /3, Dg, P e , Dg, U , and Dj. 

(3) For the initial tj and the initial iteration of the first S, estimate the starting 
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profiles of f and f with a 1/7 power law and neglect nonsimilar terms. 

(4) Calculate the coefficients a^ and h^. 

(5) Solve equations (B6) to (BIO) to obtain a new approximation for f^ and f . 

(6) Iterate steps 4 and 5 until the magnitudes of a (eq. (B2)) are approximately 
equal from successive iterations (maximum of four iterations). 

(7) Check Iq to see if [Iq| < 0. 005. 

(8) If |Iq| > 0.005, increase T) e by Arj e and go to step 4. 

(9) If |Iq| < 0. 005, the boundary- layer profiles, thicknesses, and important param- 
eters are printed out and the program proceeds to the next S (step 2). 
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Flow 



Figure 1. - Boundary-layer coordinate system along blade 
cross section. 



Ratio of boundary-layer coordinate normal to blade surface to 
boundary-layer thickness, y/6 

Figure 2. - Approximate comparison of measured mixing length with 
Prandtl's mixing length for low free-stream turbulence levels. 
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boundary-layer thickness, y/5 (b) pressure surface. 

Figure 3. - Approximate comparison of measured mixing length with Figure 4. - Comparison of predicted displacement and momentum thicknesses 

Prandtl's mixing length for high free-stream turbulence levels. with experimental data (ref. 12) along the cascade airfoil section for two- 

dimensional flow with a cascade inlet-flow angle of 55° (2D-55 0 ). 




Theory Data 
O 


Parameter 

Displacement thickness 
Momentum thickness 


Velocity gradient 
Skin-friction coefficient 


.006 1 — 
.004 — 
.002 — 
0 — 
.002 — 


.006 
.004 — 
.002 — 
0 — 
-.002 — 


8 ' 



Percent of airfoil distance to chord length, x/c 


(b) Three-dimensional flow with a cascade inlet -flow angle 
of 60° (30-60°). 


Figure 5. - Comparison of predicted displacement and momentum thicknesses 
with experimental data (ref. 12) along the suction surface of the cascade 
airfoil section for two inlet -flow conditions. 
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Transformed variable, t] 


Figure 6. - Comparison of predicted velocity profile with experimental data at 20 percent of chord for two- 
dimensional flow on the suction surface at a inlet-flow angle of 55°. 



O Data (ref. 12) 





Velocity profile, f v Velocity, f 



Figure 8. - Comparison of predicted velocity profile with experimental data at 98 percent of chord for two- 
dimensional flow on the suction surface at an inlet-flow angle of 65°. 



Figure 9. - Comparison of predicted velocity profile with experimental data at 98 percent of chord for 
three-dimensional flow on the suction surface at an inlet-flow angle of 60°. 
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